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Abstract 

We study high dimensional integration in the quantum model of 
computation. We develop quantum algorithms for integration of func- 
tions from Sobolev classes WJ([0, l] d ) and analyze their convergence 
rates. We also prové lower bounds which show that the proposed al- 
gorithms are, in many cases, optimal within the setting of quantum 
computing. This extends recent results of Novak on integration of 
functions from Hòlder classes. 



1 Introduction 

Since Shor's (1994) discovery of a polynomial factoring algorithm on a quan- 
tum computer, the question of the potential power of quantum computing 
was posed and studied for many problems of computer science. Most of these 
are of discrete type, while so far little was done for numerical problems of 
analysis. This field contains a variety of intrinsically difficult problems. One 
of tliem is high dimensional integration. 

To judge possible gains by a quantum computer, one first of all needs 
to know the complexity of the respective problem in the classical settings. 
The complexity of many bàsic numerical problems in the classical deter- 
ministic and randomized setting is well understood due to previous ef- 
forts in information-based complexity theory (see Traub, Wasilkowski, and 
Wozniakowski, 1988, Novak, 1988, and Heinrich, 1993). This theory estab- 
lished precise complexity rates by developing optimal algorithms, on one 
hand, and proving matching lower bounds, on the other. 

Based on such grounds, it is a challenging task to study these problems 
in the quantum model of computation and compare the results to the known 
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classical complexities, this way locating problems where quantum computing 
could bring essential speedups, and moreover, quantitatively assessing the 
reachable gain. 

In a series of papers, Novak and the author started to investigate this 
field. Their research dealt with summation of sequences and integration 
of functions. So Novak (2001) studies integration of functions from Hòlder 
spaces, using the algorithm of Brassard, H0yer, Mosca, and Tapp (2000) for 
approximating the mean of uniformly bounded sequences. Heinrich (2001a) 
and Heinrich and Novak (2001b) developed quantum algorithms for the 
mean of p-summable sequences and proved their optimality. Moreover, such 
an approach required a formal model of quantum computation for numerical 
problems, which was developed and studied in Heinrich (2001a). This way 
the bàsic elements of a quantum setting of information-based complexity 
theory were established. First ideas about path integration are discussed in 
Traub and Wozniakowski (2001). 

Integration of functions from Sobolev spaces is one of the bàsic numer- 
ical problems for which we know the complexity both in the classical de- 
terministic and randomized setting. In the present paper we study this 
question in the quantum setting. We develop a quantum integration algo- 
rithm by splitting the problem into levels, using a hierarchy of quadrature 
formulas, and this way reducing it to computing the mean of families of 
p-summable sequences. This enables us to apply the results of Heinrich 
(2001a) and Heinrich and Novak (2001b), and shows that the investigation 
of p-summable sequences was an important prerequisite to handle functions 
from Sobolev classes. We also prové lower bound which show the optimality 
(up to logarithmic factors) of the proposed algorithms. 

The contents of the paper is as follows. In section 2 we recali some 
notation from the quantum setting for numerical problems as developed in 
Heinrich (2001a). In section 3 we add some new results of general type 
which will be needed later on. Section 4 recalis known facts about sum- 
mation of sequences and provides some refinements of estimates. The main 
result about quantum integration of functions from Sobolev classes is stated 
and proved in section 5. The paper concludes with section 6 containing 
comments on the quantum bit model and a summary including comparisons 
to the classical deterministic and randomized setting. 

For more details on the quantum setting for numerical problems we re- 
fer to Heinrich (2001a), also to the survey by Heinrich and Novak (2001a), 
and to an introduction by Heinrich (2001b). Furthermore, for general back- 
ground on quantum computing we refer to the surveys Aharonov (1998), 
Ekert, Hayden, and Inamori (2000), Shor (2000), and to the monographs 
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Pittenger (1999), Gruska (1999), and Nielsen and Chuang (2000). 



2 Notation 

For nonempty sets D and K, we denote by T(D, K) the set of all functions 
from D to K. Let F Ç JP"(D, K) be a nonempty subset. Let K stand for 
eitl·ier R or C, the field of real or complex numbers, let G be a normed space 
over K, and let S : F — > G be a mapping. We seek to approximate S(f) for 
/ G -F by means of quantum computations. Let üi be the two-dimensional 
complex Hilbert space C 2 , with its unit vector basis {eo,ei}, let 

# m = H x <g) ••• <g)iïi, 



equippeded with the tensor Hilbert space structure. Denote 

Z[0,iV) := {0,...,N-1} 

for N G N, where we agree to write, as usual, N = {1, 2, . . . } and No = 
N U {0}. Let C m = {\i) : i G Z[0, 2 m )} be the canonical basis of H m , where 
\i) stands for ej Q <8> • ■ ■ <8> ej m _i , i = Sfc^o 1 Jfc2 m_1_fc the binary expansion of 
i. Let U(H m ) denote the set of unitary operators on H m . 
A quantum query on F is given by a tuple 

Q = (m, m' , m", Z, r, 0), 

where m,m',m" G N,m' + m" < m, Z Ç Z[0, 2" 1 ') is a nonempty subset, 
and 

r : D 
P :K^Z[0,2 m ") 

are arbitrary mappings. We let m(Q) := m be the number of qubits of Q. 

Given a query Q, we define for each / G F the unitary operator Qf by 
setting for \i) \x) \y) G C m = C m / ® C m " (g) C m _ m ^ m «: 

0,1^1^1^-/ K)k©/?(/(r«)))|y) if t€Z 
V/ W \x) \y) - | ^ |^ |^ otherwise, 

where © means addition modulo 2 m ". 

A quantum algorithm on F with no measurement is a tuple 

A=(Q,(U^ =0 ). 
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Here Q is a quantum query on F, n G No and Uj G U{H m ) (j = 0, . . . , n), 
with m = m(Q). Given / G F, we define ylj G U{H m ) as 

Af = U n Q f U n - 1 ...U 1 QfU . 

We denote by n 9 (.A) := n the number of queries and by m(A) = m = m(Q) 
the number of qubits of A. Let (^4/(x, y)) x ,j/ez[o,2 m ) be the matrix of the 
transformation Af in the canonical basis C m . 

A quantum algorithm from F to G with measurements is a tuple 

A = {{A í f í zl{b í )tlM, 

where k G N, A(_ (£ = 0, . . . , k — 1) are quantum algorithms on F with no 
measurements, 

b G Z[0,2 m °), 
6/ : ]JZ[0,2 mi ) -► Z[0,2 mí ) (1 < £ < fc - 1), 

i=0 

where := m(A^), and 

fc-i 



ip : []Z[0,2 m ') -► G. 



The output of A at input / G F will be a probability measure A(f) on G, 
defined as follows: First put 

PA,f{x ,...,x k -i) = \A j(x ,b )\ 2 \A 1 j(x 1 ,b 1 (x ))\ 2 ... 

■ ■ ■ \Ak-i,f{x k -i,b k -i(x , . . . ,x k - 2 ))\ 2 - 

Then define A(f) by setting for any subset C Ç G 

A{f){C)= ^2 pajÍxq,...^^). 

ip(xo,...,x k - 1 )ec 

Let n q (A) := Yle=o n q{Ae) denote the number of queries used by A. For 
more details and background see Heinrich (2001a). Note that we often use 
the term 'quantum algorithm' (or just 'algorithm'), meaning a quantum 
algorithm with measurement(s). 

If A is an algorithm with one measurement, the above definition simpli- 
fies essentially. Such an algorithm is given by 

A=(A ,b ,f), Ao = (Q,(Ujïï =0 ). 
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The quantum computation is carried out on m := m(Q) qubits. For / € F 
the algorithm starts in the state | òo) and produces 

\il>f) = U n Q f U n - 1 ...U 1 Q f U \bo). 

Let 

2 m -l 

(referring to the notation above, we have üíj = Aoj(i,bo)). Then A out- 
puts the element (p(i) 6 G with probability |aj,/| 2 - It is shown in Heinrich 
(2001a), Lemma 1, that for each algorithm A with k measurements there is 
an algorithm A with one measurement such that A(f) = A(f) for all f E F 
and A uses just twice the number of queries of A, that is, n q (A) = 2n q (A). 
Hence, as long as we are concerned with studying minimal query error and 
complexity (see below) up to the order, that is, up to constant factors, we 
can restrict ourselves to algorithms with one measurement. 

Let 9 > 0. For a quantum algorithm A we define the (probabilistic) error 
at / € F as follows. Let Ç be a random variable with distribution A(f). 
Then 

e(S, A, f, 6) = inf {e > | P{ \\S(f) - Q > s} < 6} 

(note that this infimum is always attained). Hence e(S,A,f,9) < e iff the 
algorithm A computes S(f) with error at most e and probability at least 
1 - d. Trivially, e(S, A, f, 6) = for 9 > 1. We put 

e(S,A,F,8) = su P e(S,AJ,e) 
fer 

(we allow the value +oo for this quantity). Furthermore, we set 
eUS,F,9) 

= inï{e(S, A, F, 9) \ A is any quantum algorithm with n q (A) < n}. 

It is customary to consider these quantities at a fixed error probability level: 
We denote 

e(S,A,f) = e(S,A, /, 1/4) 

and similarly, 

e(S,A,F) = e(S,A,F,l/4), e%(S,F) = e^(S, F, 1/4). 
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The choice = 1/4 is arbitrary - any fixed 9 < 1/2 would do. The quantity 
en(S,F) is central for our study - it is the n-th minimal query error, that 
is, the smallest error which can be reached using at most n queries. Note 
that it essentially sufnces to study ek{S, F) instead of ek{S, F, 6), since with 
0(y) repetitions, the error probability can be reduced to 2~ v (see Lemma || 
below). 

The query complexity is defined for e > by 
comp^S, F) = 

min{n ç (^4) | A is any quantum algorithm with e(S, A, F) < e} 

(we put compf (S, F) = +oo if there is no such algorithm). It is easily 
checked that these functions are inverse to each other in the following sense: 
For all n G No and e > 0, ek{S, F) < s if and only if compfj (S, F) < n for 
all Ei > e. Hence it sufnces to determine one of them. We shall principally 
choose the first one. 



3 Some General Results 

Let ^ F Ç F(D,K) and / F Ç F(D,K), where D,D,K,K are 
nonempty sets. Suppose we want to construct an algorithm A on F by 
the help of some reduction to an already known algorithm A on F in the 
following form: For / G F we construct a function / = T(f) £ F and 
apply A to it. When does this indeed give an algorithm on F? To clarify 
the problem, note that by defmition, an algorithm A on F can only use 
queries Q on F itself, while in the approach above we use Qr(f) instead, 
where Q is a query on F. The way out is to simulate Qr(f) as where B 
is an algorithm without measurement on F. The following result contains 
sufficient conditions and is a generalization of Lemma 5 of Heinrich (2001a). 

Assume that we are given a mapping r : F — > F of the following type: 
There are k, m* S N and mappings 



D^D (j = 0,...,«-l) 
K -► Z[0,2 m *) 



D x Z[0,2 m *) K -» K 
such that for f E F and s £ D 

(T(f))(s) = q(s, o f o Vo (s), ...,/3 o f o ^ 00). 
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Lemma 1. For each quantum query Q on F and each mapping T of the 
above form (^) there is a quantum algorithm without measurement B on F 
such that n q (B) = 2k and for all f G F, x G Z[0, 2 m ), 

(Qr(/) \x)) |0) m _^ = B f \x) |0) m _~ , 

where fh = m(Q), m = m{B) > fh and |0) _~ stands for the zero state of 
the last m — fh qubíts. 

Proof. Let 

Q = (m,fh',m", Z,t,P), 

and put 

kq = [log «] , m = m + kq + nm*, m! = fh' + Ko, ?n" = m*, 
Z = Zx[Q,k), r(i,j) = m(T(i)) for (i, j) G Z, 

let /? be as above, and define 

Q = (m,m' ,m" , Z,t, f3). 

We represent 

H m = H m > ® H m 'i <g> H m _ m /_ m ii (g) iï reo (8) Jï®?, 
a basis state of which will be written as 

K) k) to) li) ko) • • • • 
Define the permutation operator Pq by 

Pq \i) \ x ) \v) li) No) • • • \z K -i) = \i) li) \zq) . . ■ \z R -x) \x) \y) , 

another permutation operator 

P \i) \j) \z ) ... \zj) . . . \z K -x) \x) \y) = \i) \j) \zj) . . . \z ) . . . \z K -i) \x) \y) , 

the following counting operators 

C |») |i) . . . \y) = ...\y) 
C\i)\j)...\y) = |»>|iel)...|y), 

where © is addition modulo 2 K °, and the operator of sign inversion 
J\i)\j)\z )...\y) = \i)\j)\ez )...\y), 
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where © is subtraction modulo 2 m * and Qz stands for © z. Finally, let 
T\i) |i)|z )...|^i)|x)|y) 
= I*) li) W) ■ ■ ■ \ z k-i) x<S>f3oQ(T(i),zo,...,z K -i)^\y) 

if i G Z, and 

T \i) |i) \z ) . . . \z K ^) \x) \y) = \i) |i) \z ) . . . |z K _i) \x) \y) 
if i Z. We define B by setting for f E F, 

B f = Pu 1 Cü 1 (CQ f JP) K T(PQ f C- 1 ) K C P . 
Let us trace the action of Bf on 

|·>|s>|y>|0>|0>...|0>. 
First we assume i € Z. The application of Pq, followed by Cq, gives 

|i)|Kmod 2 K °)|0)...|0)|x)|y). 
The transformation PQfC" 1 leads to 

|i)|/e-l)|0)...| ) 3o/o7/ (S _ior(i))|a:)|y), 
and after the remaining n — 1 applications of PQfC^ 1 we get 

|«> |0) |/3 o / o rjo o f (i)) . . . |/3 o / o r/ K _i o r(i)) \x) \y) . 
Then the above is mapped by T to 

\i) |0) |/3 o / o r/o o r(«)) . . . |/3 o / o r/ K _i o r(i)) 
xffi /3o q(t(í),Po forjo °r(i), ■ ■ ■ ,/3° / o 77 K _i or(i))^ |y) 
= K) 1°) 1/3 ° / ° % ° ^(0) • • • \P ° / ° Vk-i ° 

xe)8((r(/))(r(t))))|y>. 

The transformation (CQfJP) K produces 

|»> |« mod 2 K °) |0) . . . |0) |x © £((r(/))(r(i)))) |y) , 
and finally Pq -1 ^ 1 gives 

|i) |x © j9((T(/))(r(i)))) |y> |0) |0> . . . |0) = (Q r(/) \i) |x) \y)) |0) m _~ . 
The case i Z is checked analogously. □ 



Corollary 1. Given a mapping F : F — » F as in a normed space G 
and a quantum algorithm A from F to G, there is a quantum algorithm A 
from F to G with 

n q (A) = 2nn q {A) 

and for all f G F 

A(f) = A(T(f)). 

Consequently, if S : F — > G is any mapping and S = S o F, then for each 
n G No 

e q 2Kn (S,F)<el(S,F). 

The proof is literally the same as that of Corollary 1 in Heinrich (2001a). 
We omit it here. 

Lemma 2. Let D,K and F Ç !F(D,K) be nonempty sets, let k G No 
and let Si : F — ► R (l = 0, . . . , k) be mappings. Define S : F — > R 
#(/) = Ef=o^(/) (f ^ F )- Let °o, ■■■ A > 0, n , ... G N and put 
n = Y2i=o n i- Then 

k k 

el{S,F^i)<Y. e n^ F ^)- 

1=0 1=0 

Proof. Let 5 > and let A\ be a quantum algorithm from F to R with 
" g (^z) < n\ and 

Let A = J^f = o -A/ be the composed algorithm (in the sense of section 2 of 
Heinrich, 2001a). Then 

k k 

n q (A)=Y j n q (A l )<^2n l . (2) 

1=0 1=0 

Fix an / G F and let (0,/)f=o be independent random variables with distri- 
bution Ai(f). It follows that with probability at least 1 — 0j, 

|5K/)-Cd<4 ! (^^) + <5· 

Setting 

fe 

C/ =2^G,/i 
z=o 



9 



we infer from Lemma 2 of Heinrich (2001a) that C/ has distribution A(f). 
Consequently, 



\s(f) - C/1 = | " ^ EW 5 <' F ^) + *) 

with probability at least 



1=0 1=0 



z=o 

This gives 

fe fc 
z=o z=o 

for all / 6 .F, and the desired result follows from (||). □ 

Corollary 2. Let D,K, F C F{D,K), k £ N and S, S", : F -> R (/ = 
0, . . . , k) be as in Lemma [|. Assume uq, . . . , E N satisfy 



k 1 

y e -^/ 8 <i. 

- 4 

z=o 

Leí uq, . . . ,nk £ No an<i piií n = Xw=o ^í 7 "^- 

4(^)^4,(5^). 

1=0 

This is an obvious consequence of Lemma |2| above and of Lemma 3 in 
Heinrich (2001a), which can be restated in the following form: 

Lemma 3. Let S be any mappíng from F Ç J-(D, K) to R. Then for each 
n, v G N, 

el n (S,F,e-»/ 8 )<el(S,F). 
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4 Summation 



This section provides the prerequisites from summation needed later for the 
study of integration. For N € N and 1 < p < oo, let Lp denote the space 
of all functions / : Z[0, N) — > R, equipped with the norm 

/ jv-i \ Vp 

11/11^ = [n£ i/«i p J • 

Define 5at : Lp —> R by 

1 

and let 

S(^) :={/€L^|||/||^<l} 

be the unit ball of Lp . 

We need the following results about summation, where @ and @ are 
from Heinrich (2001a), Theorems 1 and 2, and (||) is from Heinrich and 
Novak (2001b), Corollary 2. 

Proposition 1. Let 1 < p < oo. There are constants 01,02,03 > such 
that for all n,N G N wí/í n < ciJV, 

C2H- 1 <el(S N ,B{L^)) <c 3 n~ 1 if 2 < p < 00, (3) 
c 2 n~ 1 < el(S N ,B(Lf)) < c^ 1 log 3 / 2 nloglogn, (4) 

and 

c 2 min(n- 2 ( 1 - 1 / p ),n- 2 / p iV 2 /P- 1 ) < e q n {S N ,B(L%)) 

< c 3 min(n- 2(1 - 1/p) ,n- 2/p iV 2/p - 1 )max(log(n/VÏV),l) 2/p " 1 . (5) 

*/ 1 < p < 2. 

Remark. We often use the same symbol c, ci, . . . for possibly different pos- 
itive constants (also when they appear in a sequence of relations). These 
constants are either absolute or may depend only on p,r,d - in all lemmas 
and the theorem this is precisely described anyway by the order of the quan- 
tifiers. 
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In the case p = 2 we will not use Proposition || alone - that would give 
just a logarithmic factor instead of the iterated logarithm of Theorem [ï] be- 
low. In the region where n is close to N we use a refinement which can 
be obtained on the basis of the results in Heinrich and Novak (2001b). We 
introduce for M £ N 

SnmÍ = jf Yl /(*) 
ieZ[0,N], \f(ï)\<M 

and 

S N,Mf = S NÍ ~ S N ,MÍ = jx f®' 

ieZ[0,N], |/(í)|>Af 

Let us first recali the case p = 2 of Corollary 3 of Heinrich and Novak 
(2001b), which we will use here: 

Lemma 4. There is a constant c > such that for all n, M, JVéN, 
whenever 

M > cNn' 1 max(log(n/v / ÏV), 1). 

The next result can be shown by repeating the respective part of the 
proof of Theorem 1 of Heinrich (2001a) (for the analogous p < 2 case see 
also the proof of Proposition 2 in Heinrich and Novak, 2001b). 

Lemma 5. There is a constant c > such that for all k,n, N £ N , k > 1 
(%(S Nt2 k,B(L?)) <c(n- 1 k 3 / 2 logk + 2 k n~ 2 (klogk) 2 ^j . 

From these we can derive the following estimate: 

Lemma 6. There is a constant c > such that for all n, N £ N with 
n< N, 

e« n (S N ,B(L%)) < cn" 1 A(n,iV) 3 / 2 logA(n,iV), 

where 

A(n, N) = log(N/n) + log log(n + 1) + 2. 

Remark. Observe that Lemma || gives an improvement over Proposition [ï] 
only for n close to N. 
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Proof. It suffices to prové the statement for 

N < ra 3/2 , (6) 

the other case follows directly from (Q). Let co denote the constant from 
Lemma |] and let fe be the smallest natural number with fe > 2 and 

CoNn- 1 max(log(n/VÏV), 1) < 2 k . (7) 

By Lemma |], 

el(S> N2k ,B(L»)) = 0. (8) 
This together with Lemma || and Corollary |2| gives 

e q Cin (S N ,B(L$)) < cn- 1 ^ 2 logk + c2 k n- 2 (klogk) 2 , (9) 
with a certain constant c\ £ N. It follows from (^) that 

2 k ~ x < max (coNn- 1 max(log(n/\/ÏV), 1), 2), (10) 
which, in turn, implies 

2 k < cNn- 1 log(n+ 1), (11) 
fe < c(log(iV/n) +loglog(n+ 1) + 1) = c\(n,N), (12) 

and thus 

k < clog(N + 1). (13) 
From (U) and A(n, N) > 2 we conclude 

logfe < clogÀ(n, iV), (14) 

while ([[|) gives 

log fe < c(loglog(iV + l) + 1). (15) 

From (0), (||), (H), and (§), we infer 

2 fc n _1 fe 1/2 (logfe) 2 < cAfn" 2 log(n + l)(logAf) 1/2 (loglog(iV + l) + l) 2 
< cA^ 1/3 (log(iV + l)) 3/2 (loglog(iV + l) + l) 2 <c. 

Consequently 

2 fc n~ 2 (felogfe) 2 < cn^k 3 / 2 logk, 
and hence, by (|), ©, and (0), 



e2 in (5jv,jB(L^)) < cn^fe 3 / 2 log fe < cn _1 A(n, iV) 3/2 log A(n, iV). 



□ 
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5 Integration 



This section contains the main result. Let D = [0, \] d and let C(D) be the 
space of continuous functions on D, equipped with the supremum norm. For 
1 < p < oo, let L p (D) be the space of p-integrable functions, endowed with 
the usual norm 

Wf\\L p (D ) =( K l D \f(tWdt 

The Sobolev space W£ (D) consists of all functions / € L p (D) such that for 

all a = (ai, . . . , a^) £ w ^ tn M := Sj=i a i — r > the generalized partial 
derivative d a f belongs to L p (D). The norm on W£ (D) is defined as 

Wf\\w;(D) = i E W^fWliD) 

\\a\<r 

We shall assume that r/d> l/p, which, by the Sobolev embedding theorem 
(see Adams, 1975, or Triebel, 1995), implies that functions from Wp(D) are 
continuous on D, and hence function vàlues are well defined. Let B(Wp(D)) 
be the unit ball of W p (D) and let l d ■ Wp(D) — ► R be the integration 
operator 

W) = i f(t)dt. 
Jd 

Theorem 1. Let r, de N, 1 < p < oo and assume r/d > l/p. There are 
constants c±, c<i > such that for all n £ N with n > 4 

cm-^ 1 < e£(/ d ,£(W£ (£>))) <c 2 rr r / d - 1 i/ 2 < p < oo, 
Cl n- r ' d ~ l < el(I d ,B(W^D)))<c 2 n- r / d - 1 X (n), 

Cin -r/d-i < e q n (I d ,B(W;(D))) < c 2 n- r / d - 1 {\ogn) 2 l p - 1 if 1 < p < 2. 

T/ie function Ào denotes an iterated-logarithmic factor: 

Xo(n) = (loglogn) 3 / 2 logloglogn. 

Proof. First we prepare the needed tools for the upper bound proof. For 
l G N let 

£>= u D « 
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be the partition of D into 2 dl congruent cubes of disjoint interior. Let su 
denote the point in Du with the smallest Euclidean norm. We introduce the 
following extension operator 

E U :F(D,-R) ^F(D,K) 

by setting 

(E li f)(s) = f(s H + 2- l s) 
for / G F{D, R) and s G D. Now let J be any quadrature rule on C(D), 

K-l 

3=0 

with a,j G R and tj G D, which is exact on V r -\{D), that is, 

J/ = / d / for all f£V r -i{D), (16) 

where V r -i(D) denotes the space of polynomials on D of degree not exceed- 
ing r — 1. (For example, for d = 1 one can take the Newton-Cotes formulas 
of appropriate degree and for d > 1 their tensor products.) Since r > d/p, 
we have, by the Sobolev embedding theorem (see Adams, 1975, or Triebel, 
1995), Wp(D) C C(D) and there is a constant c > such that for each 
/ e W£(D) 

||/||C(D)<C||/||^(D)- (17) 

Consequently, 

K — 1 K — 1 

|J/| < ^KH/feoi <EKIH/IIc(d) < cII/Hh-p). (is) 

3=0 j=0 

For / G Wp (D) we denote 

. Vp 

I/|,, p ,d= í £ J D \d a f(tWdt 

\\a\=r D j 

According to Theorem 3.1.1 in Ciarlet (1978), there is a constant c > such 
that for all / G W;{D) 

H-f ~ 9\\w;(D) < c\f\ r ,p,D- (19) 
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We conclude from (0), fljj) and @, 

IW-J/I < inf \I d (f- g )- J(f-g)\ 
geV r -i(D) 

< c inf \\f-g\\ W r( D )<c\f\ r , p , D . (20) 

geV r - 1 (D) 

Now define for / G No 

2 dl -l 2 m -1k-1 

J,/ = 2- dí ^ J(E H f) = 2- dl £ ^a J /( S/ï + 2-%), 

i=0 j=0 j'=0 

which is the composed quadrature obtained by scaling J to the subcubes 
Du. Then we have for / G 

2 di -l 



IW- Jí/| = |W-2- dí ^ J{E H f)\ 

i=0 

2 m_ 1 

< 2- dl \I d (E u f)-J(E u f)\ 

í=0 
2 di -l 

< c2~ dl Y, \Elif\r, P ,D 

i=0 

I 2 di -l 

< c r*J]|^ 



i=0 



and 

2 tí(_ 1 2 dl -l 



i=0 ï=0 \a\=r JD 

2 dl — 1 

2-pri í i^/wr^ 



i=0 |a|=r" 

= 2-^|/l^ < 2- pri ||/||^ (D) . (21) 

It follows that 

|W - Jif\ < c2- rí |/lr, P ,D < c2- w ||/||^ (Z3) . (22) 
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Let us now describe the main idea: First we approximate Idf by the 
quadrature J^f for some k, giving the desired precision, but having a num- 
ber of nodes much larger than n. This in turn, will be split into the 
sum of a single quadrature J^ , with number of nodes of the order n, which 
we compute classically, and a hierarchy of quadratures (more precisely, dif- 
ferences of quadratures) J[ (l = kç,, . . . , k — 1). It will be shown that the 
computation of the J[f reduces to the computation of the mean of sequences 
with well-bounded L^-norms for suitable N\. This enables us to apply the 
results of section 4 and approximate the means by quantum algorithms. In 
the sequel we give the formal details, the proper balancing of parameters 
and the proof of the error estimates. 

Define 

2 d -lK-l K-l 

J'f:=(Ji-J )f = 2- d ^^a j /(s M + 2- 1 í j )-^a J /(í j ) 

i=0 j=0 j=0 

k'—1 

■= E4/(^)' ( 23 ) 

3=0 



where 



For Z G N , set 



k' < n(2 d + 1). (24) 



k'-I 

J' u f = j'(E li f)=Y,a j f(su + l~%), (25) 

3=0 

2 dl -l 

j> = 2 -dl ^ j>.. (26) 

i=0 

It is easily checked that 

2 dl -l 

J l+1 f = 2~ dl Ji(Euf). 

i=0 

and hence 

2 dl -l 

Ji+if-Jif = 2- dl Y,(-h(Euf)-Jo{Euf)) 

i=0 
2 dl -l 

= 2- dl J'uf = Ai- (27) 

i=0 
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Using (gg) and fl2T|), we get 

2 dí_ 1 2 di -l 



2^ £ \W ^ 2- dl Y,\UEuf)-MEuf)\ P 

i=0 i=0 
2 dl -l 

< 2- dl E (K/d - Ji)(E H f)\ + \{i d - J )(E h )fW 



i=0 
2 dl -l 

< c2- dl £ < c2-^||/||^ r(D) . (28) 

Now we derive the upper bounds. Clearly, it sufíices to prové them for 

n > max(ít, 5). (29) 

Let 

k =[d- 1 log(n/ K )\. (30) 
By the above, we have k$ > 0. Furthermore, let 

k=\{l + d/r)ko], (31) 

hence k > k . By @) 

k-l 

Jk = Jko + J2 J í- ( 32 ) 
l=k 

For 

k <l <k (33) 

put iV; = 2 dí . We shall define mappmgs T/ : B{W r p {D)) -» in order to 
apply Lemma |]. For this purpose we fix an m* £ N with 

2 -m*/2 < k -l 2 ~ rk (34) 

and 

where c is the constant from ([Ï7|). Hence, 

\ C(D) < 2 m */2-i for /eB (W;p)). (36) 
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Define Vlj : Z[0, JV,) D (i = 0, - 1) by 

%·(*) = «íi + 2- í ^ (i€Z[0,JV,)), 

and : R Z[0, 2 m * ) for z 6 R by 

r if * < -2 m */2-i 

= i [2 m * /2 (z + 2™*/ 2 " 1 )] if - 2 m *l 2 - 1 <z< 2 m */2-i 
[ 2 m * - 1 if z > 2" 1 */ 2 - 1 . 

Furthermore, let 7 : Z[0,2 m *) -> R be defined for y G Z[0,2 m *) as 

7 (y) = 2 " m */2 y _ 2 m*/2-l_ 

It follows that for _2 m */2-i < z < 2^/2-^ 

7(/3(^)) <^<7(/3W)+2- m * /2 . 
Next let q : Z[0,2 m *) K ' -> R be given by 

Q(yo, Vk'-x) = a í^(%')- 
3=0 



Finally, we set 

r i (/)« = í?((/?o/°%·«)·=o 1 )· 

for / G B(W£(.D)). We have 

r z (/)(z) = ^4 7 W(^ + 2-^))), 
3=0 

hence, by fl2|), (||) and (||), 

|Jíí/-W)(*)| < ^|a;.||/( s/ï + 2- í í;.)-7(/3(/(^ + 2-'í;-; 



3=0 

k'-1 

< 2- m ^ 2 ^2\a' 3 \ <c2- m *' 2 Kck^T 
3=0 



rk 



and therefore, by @, for all / G B(W£(D)), 

2 dl -l 

\j[f - s^vm < 2- di Y,\ j u- wxoi < ck-'2 



i=0 
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Moreover, by (I 



and (p 



\m)\\ L ^ < \\(Jíif)^a\\ L »i + \\W)-(Jíif) 



i=0 W L N i 



< 11(4/) 



l f sN t -l 



i=0 11^ + IIW) - {J'uf) 



i=0 W L N i 



c2 



Consequently, 

T l (B(W;(D)))Çc2- rl B(L^). 
By @), K2 dko < n, hence 



(42) 



(43) 



(this just means that with n2 dk ° queries we can compute Jfc , the mean 
of n2 dk ° numbers, classically, or, more precisely, up to any precision by 
simulating the classical computation on a suitable number of qubits). By 
assumption, r/d > 1/p and p > 1. Hence 



d 

r > - > 
P 



P 



1 d. 



Now fix any 5 with 



< d < min ( r, — ( r 



- - 1 ] d 

P 



(44) 



and put for l = kç, , . . . , k — 1 



ni 



2<ifco— 5(1— kg) 



(45) 



vi = [8(2 ln(Z - fc + 1) + In 8)] . (46) 
It follows from fl46]) that 

^^<^(l-fc„ + l)- 2 <i (47) 
i=feo i=fco 

Put 

jfc-1 

n = n + 2/í' z^ïfy. (48) 

l=k 
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By <m, <m, and (B 



k-l 



n 



< u + 2k'Y^ [8(2 ln(Z - fc + 1) + hi 8)] 2 d/c °- <5 ( / - fc °) 
í=feo 



< c2 dfc °<cn. 



From above and Lemma 6 (i) of Heinrich (2001a), 



el(I d ,B(W;(D)))<c2- rk + el(J k ,B(W;(D))). 



(49) 



(50) 



By Lemma y and 

4(J fc> iB(WÇ(D))) 



< eUJk ,B(W;(D)),0) + eí_ n (J k - J hv B{W;{D))) 

= 4- n ( J k-Jk ,B(w;(D))). 



(51) 



Using (H), (pi), Corollary | and 



4_JJ fc - J fco ,fí(^;(£))) 



we get 



fc-i 



^J/,B(WJ(D)) 



J=feo 



fc-1 



Z=fc 



(52) 



From ([ïï|) above and Lemma 6 (i) of Heinrich (2001a) we conclude 

eLn ! (^^(% r P)))<^- 1 2^ fc + eL ni (%r^(WÇp))). (53) 
Corollary [ï], relation fl42] ) above and Lemma 6(iii) of Heinrich (2001a) give 



e« (S^S^I)))) < e* ni (S Nl ,c2- rí B(Lp)) 



c2- rl eUS Nl ,B(L^)). (54) 



Joining (p0|)-(p4^, we conclude 



fc-1 



e|(J d , B(WJ (2>))) < c 2-^ + c £ 2"'^ , B(L?)). (55) 

Z=fc 
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Now we prové the upper bound in the case 2 < p < oo. Relation 
Proposition @, ©, and © 



k-l 



-W -1 



n, 



l=k 



k-l 



< c2~ rk + c 2~( r+d ï k ° ^2 2~( r ~ s ï( 1 ~ k °} 

l=ko 

< c 2-( r+d ^ < cnr r l d - 1 . 



Next we consider the case 1 < p < 2. Observe that, by (44) 



-6 < r - (- - 1 ] d. 

p \p 



It follows from (H), Proposition 0, f)§), Q), and @ that 



4(I d ,B(W;(D))) 



k-l 



< 



c2- rk + cY,l~ rl n7 2/P rf /p - 1 



max 



(log^/v^),!) 2 /^ 1 



l=ko 
k-l 



< c2~ rk + c 2- rl -í dk0+ í 5il - ka)H í- 1)d \k + l)^ 1 



l=k 



" 1 ^ (-r+(|-l)«H-|i)(i-fco) 



< c2- rfc + c2-( r+ ^ fco (A;o + l)i ; " 1 J^2 

z=fc 

< c2-( r+d ) fco (A;o + < cn~ r/,i ~ 1 (logn)Í~ 1 . 

Finally we consider the case p = 2. From (|55|) and Lemma || we get 

eí(Id,B(W;(D))) 

k-i 

< c2~ rk + c ^ 2- ri n- 1 A(n z , iVO 3 / 2 log \(n h N t ) 



l=ko 

< c2- rk + c2' i - r+d ^ ko x 
k-i 

2-ir-W-ko) (| _ A; + log log n + 1) 3/2 log(/ - A; + log log n - 

í=feo 

< c2~( r+d ) feo (loglogn) 3 / 2 logloglogn < crT r/d ' 1 \ {n). 
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To conclude the proof of the upper bounds in all three cases, we use (49) 
and scale n. 

Now we turn to the lower bounds. Since B(W£(D)) C B(W^(D)) for 
p > q, it suffices to consider the case 2 < p < oo. Fix such a p. Let ip be a 
C°° function on ü d with 

supp^ C (0,l) d , ai:=/ d V>0, 

and denote ||V ; llvK r (D) = a 2- Let n € N, A; = [~cf _1 (log(n/ci) + l)] , where c\ is 
the constant from Proposition [T], which can be assumed to satisfy < c\ < 1, 
and put N = 2 dk . It follows that 

Ci2 -d 2 dfc < 2n < ci2 dfc = ciiV. (57) 

Set 

^(í) = V(2 fc (í-s i )) (i = 0,... J JV-l), 
with the Sj as in the beginning of the proof. We have 

I^i = 2- dk I d iJ; = a x 2- dk = a x N~ l (58) 

and 

Ui\\wr(D) < 2^5*11^11^) = a 2 2^ d M k . 

Consequently, taking into account the disjointness of the supports of the ipi, 
for all cu G R (i = 0, . . . , N - 1), 

N-l N-l 

II E a ^\L = E i«*ni^il^( D) < ^"líoOSIV (59) 

i=0 "'p^ i=0 

Fix any m* E N with 

m*/2 - 1 > dk/p. (60) 



Let /? : R Z[0, 2 m * ) and 7 : Z[0, 2 m * ) — > R be defined as in @ and 
For / G B{Lp) we have 

|/(*)l < N 1/p = 2 dfe/P < 2" 1 */ 2 " 1 . 

Hence, by (||), 

tW(0)) < /(*) < 7(/W))) + 2" m * /2 . (61) 
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Define 



N-l 



T:B{L%)^W;{D) by r(/) = £ 7 ° /? o f(i) ^. 



t=0 



By © and ©, for / e B(L^), 

||r(/)|| W;(D) < o/Jo/IU? 

< ^2 ï · fc (||/||^ + ||/- 7 o/3o/||^ 

< a 2 2 rk (l + 2 

Furthermore, by (|58|), 



AT-l 



i=0 

AT-l 

= axS^d o /3 o /). 

Define 

77 : D -» Z[0,iV) by ??(s) = min{i | € D ki }, 
with the Dfcj as in the beginning of the proof, and 

g:í}xZ[0,2 m VR by g(s, z) = j{z)^ v{s) (a). 

Then 

= l°P°f(v(*))il>n(8)(8) 
= q(s,(3o forj(s)). 

So r is of the form (^|) (with k = 1) and maps 

B(L%) into <j 2 2 rk (l + 2- m *l 2 )B{W;(D)). 
By Corollary [ï] and Lemma 6(iii) in Heinrich (2001a), 

eUh oT, B{I$)) < é n (l d ,a 2 2 rk (l + 2- m */ 2 )B(W;(D)) 
= a 2 2 rk (l + 2- m *l 2 )el{I d ,B{W;{D))) 
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Using (pH) again, we infer 



sup \S N f-S N { 1 o(3of)\<2- m */\ 

and hence, by Proposition [ï] and Lemma 6(i) and (ii) of Heinrich (2001a), 
using also (62) and (|57|), 



cn 



< el n (S N ,B(L»)) 

< eUS N o 1 op,B{L^)) + 2- m */ 2 
= aï l el n (I d or,B(L»)) + 2- m */* 



-m*/2 



el{I d ,B(W;{D))) + 2- 



-m*/2 



< cn r / d eï l (I d ,B(W;(D))) + 2- m */ 2 . 
Since m* can be made arbitrarily large, the desired result follows. 



□ 



6 Comments 

Let us discuss the cost of the presented algorithm in the bit model of com- 
putation. The algorithm consists of quantum summations on the levels 
l = ko,...,k- 1. On level l we have N t = 2 dl and n x = O^-^-^n), 
where S > does not depend on l or n. Recali also that 2 dk ° = B(n) and 
k — k() = 0(logn). Referring to the respective discussion of summation in the 
bit model in Heinrich (2001a), section 6, and Heinrich and Novak (2001b), 
section 5, we conclude that on level l we need C(log N{) qubits, 0{n\ logiVj) 
quantum gates, C(logn/ loglogn/) measurements in the case p > 2 and 

OCnfJVfVmaxílogCni/V^), 1) + log(AV^) loglog(JV,/n,)) 

measurements for p < 2. Summarizing, we see that altogether the algorithm 
needs O(logn) qubits, O(ralogn) quantum gates, C((logn) 2 loglogn) mea- 
surements for p > 2 and C(n/logn) measurements for p < 2. Thus the 
quantum bit cost differs by at most a logarithmic factor from the quantum 
query cost 0(n). 

In the following table we summarize the results of this paper and com- 
pare them with the respective known quantities of the classical deterministic 
and randomized setting. We refer to Heinrich and Novak (2001a) and the 
bliography therein for more information on the classical setting. The respec- 
tive entries of the table give the minimal error, constants and logarithmic 
factors are suppressed. 
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deterministic 


random 


quantum 


B(W; 4 (D)),2<p<oc 


n -r/d 


n -r/d-l/2 


n -r/d-l 


B(W; 4 (D)),l<p<2 


n -r/d 


n —r/d-l+l/p 


„— r/d— 1 


B{W[ d {D)) 


n -r/d 


n -r/d 


„— r/d— 1 



The quantum rate for 1 < p < 2 is a certain surprise. Previous results led 
one to conjecture that the quantum setting could reduce the exponent of the 
classical randomized setting by at most 1/2. Now we see that in the case 
p = 1 there is even a reduction by 1. 
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